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ABSTRACT 


v 


This report covers work done under NASA Grant NGR 09-015-002 Supplement 
No. 57. An optical transfer function has been computed for the retroreflector array 
carried by the Starlette satellite (1975 10A). The range correction is given for extrap- 
olating laser range measurements to the center of mass of the satellite. The gain 
function and active reflecting area of the array have been computed for estimating 
laser-echo signal strengths. 


v 



OPTICAL TRANSFER FUNCTION 

OF STARLETTE RETROREFLECTOR ARRAY 

Technical Report 
RTOP 161-05-02 

lo INTRODUCTION 

This is the third in a series of reports giving optical transfer functions for satel- 
lites with retroreflector arrays- The first, performed under NASA Grant NGR 
09-015-196, presented results calculated for the BE-B, BE-C, Geos 1, Geos 2, DIG, 
DID, and Peole satellites. The second report gave results for the NTS-I array 
(1974 056A); that work was supported by Grant NGR 09-015-002 Supplement No. 57. 
The development of computer programs used in these analyses was begun under Grant 
NGR 09-015-164, "Use of a Passive Stable Satellite for Earth Physics Applications." 

Data on the Starlette satellite were obtained from the Centre National d' Etudes 
Spatiales, Toulouse, France. 
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2. ARRAY COORDINATE SYSTEM 

The position and orientation of each cube corner in the array are given by the six 
numbers, x, y, z, 0, <j>, a. The origin of the x-y-z coordinate system is the center of 
the satellite. The angles 0 and <j> are given in an x'-y'-z' coordinate, which is parallel 
to the x-y-z system (Figure la). The angle a is shown in Figure lb. The (3 and y axes 
point east and north, respectively; in other words, y is in the direction of decreasing 
4>, and (3 is in the direction of increasing 9. 

501-065 



a) b) 

Figure 1. Coordinate system for cube-comer orientation. 
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3. RETROREFLECTOR SPECIFICATIONS 

The cube corners used on Starlette have a circular face 32.8 mm in diameter, 
with a 23.3-mm length from the vertex to the face. The material is fused silica, and the 
reflecting faces are silver coated. The dihedral angles between the back faces are off- 
set from 90° to provide the necessary beam divergence to account for velocity aberration. 
The specifications on the angular offset of the reflected beam from the cube corners 
are as follows: 

Average offset = 7 arcsec, 

50% between 5 and 9 arcsec, 

100% between 0 and 10 arcsec. 

The beam divergence is related to the dihedral-angle offset by the formula 
0 = 4V 2/3 n& , 

where 


9 = angular offset of reflected beam, 
n = index of refraction of fused silica (1. 457), 

5 = offset of dihedral angles from 90°. 

The dihedral-angle offsets corresponding to the offsets of the reflected beam listed 
in the specifications have been calculated and are listed in Table 1. 
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Table 1. Cube-corner specification conversion. 


Offset of 
reflected beam 
(arc sec) 

Dihedral-angle 

offset 

(arcsec) 

0° 

0.00 

5 

1.05 

7 

1.47 

9 

1.89 


10 


2.10 



4, GEOMETRY OF THE ARRAY 


The Starlette retroreflector array consists of 60 cube corners mounted on a 24-cm 
ball. The core of the satellite is an icosahedron. Each of the 20 triangular faces of 
the core is covered by a spherical cap containing three cube corners (see Figure 2). 

A detailed geometrical analysis of the core is given in the document ’’Project 
STARLET - Etude Geometrique de la Structure Porteuse de Reflecteurs Laser, ” 

Centre Nationale D’Etudes Spatiales Publ. No. 645, Toulouse, France, 27 October 

1972 ' 501-065 



The cube corners are recessed below the surface by 0. 5 ± 0. 2 mm. The recess, 
together with the fact that the surface of the cube corner is flat, places the center of 
the front face of the cube corner closer to the center of the satellite than the 12-cm 
radius of the satellite. In the diagram below, 

R = satellite radius (120 mm), 

r = cube-corner radius, 

= 32. 8/2 = 16. 4 mm, 
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D = distance from, center of satellite to center of front face of cube comer 



Table 2 lists the coordinates, in meters, of the center of the front face and the 
orientation angles, in degrees, for each cube comer. The first column gives the index 
identifying each of the 20 faces, and the second indexes the three retroreflectors on 
each face. 
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Table 2. Cube-corner positions and orientations 
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5. REFLECTIVITY HISTOGRAM OF SPHERE 

The cube corners contributing to the reflected signal are contained in a cap whose 
half-angle is the cutoff angle of the cube corners. The earliest possible reflection 
would come from a reflector directly facing the incident beam, and the latest, from a 
cube corner near the cutoff angle. The apparent spread in range is the difference 
between the apparent reflection points for these two cases along the direction of illumi- 
nation. For a reflector facing the incident beam, the apparent reflection point is a 
distance nL in back of the front face, where n is the index of refraction (1. 457) and L 
is the length of the cube corner (23. 3 mm). The earliest reflection point is 

D - nL = 118. 37 - 1. 457 X 23. 3 = 84. 47 mm 

from the center of the satellite. The cutoff angle 6 of the cube corners is 

c 



For a cube comer at the cutoff angle, the distance of the center of the front face of the 
reflector from the center of the satellite along the direction of illumination is given by 

D cos 9 = 118. 37 cos 57? 27 = 64. 00 mm. 

The correction for optical path length in the reflector is 

L Vn 2 - sin 2 0 O = 23. 3 Vl. 457 2 - sin 2 57. 27 = 27. 72 mm . 

V 

Therefore, the latest apparent reflection point is 64 - 27. 72 = 36. 28 mm from the 
satellite's center. The range difference between the earliest and the latest reflection 
points is 84.47 - 36.28 = 48. 19 mm. Figure 3 shows the reflectivity of each 1-cm 
interval, in units of the equivalent number of cube comers at normal incidence, start- 
ing from the earliest reflection point. This histogram is for a direction of illumination 
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given by 0 ^ -13° and cj> = 60°, where 0 and (j> were defined in Figure 1. For this direc- 
tion, the total effective reflecting area is 3. 19587 times the area of one cube comer, 
and the mean apparent reflection point is 73. 4 mm from the center of the satellite. 

The cube comer shown on the histogram is in the closest possible position to the 
observer. (None of the cube comers points exactly at the observer for this particular 
direction of illumination. The earliest and latest apparent reflection points for this 
case are 82 . 2 and 41, 6 mm, respectively. ) Table 3 lists the percentage of the total 
return in each interval and the cumulative percentage, starting from the earliest 
reflection point. 


Table 3. Percentage of return in each 1-cm interval. 


Interval 

Percentage 
of return 

Cumulative 

percentage 

1 

44 

44 

2 

38 

82 

3 

15 

97 

4 

2 

99 

5 

1 

100 
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REFLECTIVITY (EQUIVALENT NUMBER OF CUBE CORNERS) 
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DISTANCE FROM CENTER OF SATELLITE (cm) 


Figure 3. Reflectivity histogram of Starlette. 
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6. SIGNAL-STRENGTH COMPUTATION 


The data contained in the tables presented later can be used to estimate signal 
strengths for laser ranging, by use of the following equation; 

2 

N = hv G T A S G S A R^4 71 ’ 

where 



number of photoelectrons, 
transmitted energy, 

Planck’s constant, 
frequency of laser light, 

"gain" of transmitter, 

active reflecting area of satellite, 

’'gain" of satellite array, 
area of receiving telescope, 
atmospheric-transmission factor, 
range from station to satellite, 

constant, which includes the quantum efficiency of the photomultiplier 
and the optical transmission factors of the transmitter, the satellite, 
and the receiver. 


If the transmitted beam is a uniform spot of solid angle S2,p the "gain" function of 
the transmitter is 


= _L 

T ' 
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The author is indebted to P. O. Minott (Goddard Space Flight Center) for pointing 

out that the gain functions used in this equation do not contain the factor of 4ir used in 

the standard definition of gain. Those given in later sections can be converted to the 

standard definition by multiplying by 4v. Furthermore, the signal -strength equation 

2 

above can be converted to the standard definition of gain by adding the factor 1/ (4 -it) . 
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7. METHOD OF COMPUTING TRANSFER FUNCTION 


In computing the transfer function of the Starlette retroreflector array, the cube 
corners have been modeled as isothermal, geometrically perfect reflectors (except 
for the dihedral-angle offset) with perfect metal reflecting coatings on the back faces. 
The primary effect of real metal faces is a decrease in the intensity of the return 
signal because of the triple metallic reflection. This loss should be added to the 
constant ij of the previous section, along with reflection losses at the front face as 
the beam enters and leaves the cube comers. 

The computation of the range correction includes a correction for the optical path 
length of the ray within the cube corner. The range correction is the difference 
between the centroid of the actual return signal and the centroid of the return signal 
that would be received from a point reflector at the center of mass of the satellite. 

The correction listed is the one-way correction. 

The gain functions and range corrections presented are for the incoherent case; 
that is, the intensities of the reflections are added without taking into account coherent 
interference among the reflected signals from the individual cube corners. 

The variation of the range correction due to optical coherence has been derived 
by statistical analysis of a set of coherent returns constructed by assigning random 
phases to the reflection from each cube corner by using a pseudo random-number 
generator. Because the coherent variations are computed from a limited number of 
cases, the results should be considered as only an indication of the magnitude of the 
effect. The incoherent range correction is the mean value of the coherent corrections. 
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8. ACTIVE REFLECTING AREA AND GAIN FUNCTION 


In order to calculate reflected signal strengths from the equation in Section 6, we 
need to find the active reflecting area A g and the gain function G g . The average A g 
computed from a set of 25 different directions of illumination is 3. 197 times the area 
of one cube corner at normal incidence. The root-mean- square (ms) deviation for 
the set is 0.056 times the area of one reflector. 

Table 4 gives the gain function of the Starlette retroreflector array for two wave- 
lengths and various dihedral -angle offsets. As stated in Section 3, the average offset 
of the reflected beam is 7 arcsec, which corresponds to a dihedral -angle offset of 
about 1. 5 arcs ec. The tables for 1.5 arcsec should yield the best estimate of the 
behavior of the array. The gain function for other offsets is given for the purpose of 
estimating the extremes of behavior possible within the manufacturing tolerances used 
for the cube corners . No attempt has been made to study variations in performance 
due to factors other than dihedral-angle offset. 

The first column in Table 4 is the velocity aberration, in microradians, and the 

7 

second gives the gain function of the array, in units of 10 . To the right of the columns 
is a computer-plotted graph, in which the gain function increases to the right and the 
velocity aberration increases down the page. The gain function is the average value 
around a circle in the far field with a radius equal to the velocity aberration. The 
second part of each table gives the rms variation of the gain function around the circle. 
The computations were done for a direction of illumination given by 0 - -13° and 
4> = 60° (see Figure la). We chose this direction because the active reflecting area 
and ran g e correction for this particular case are close to the average values. 
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Table 4 


’’Gain 1 ' versus velocity aberration,, (The gain functions presented here 
differ by a factor of 4 tt from the standard definition of gain; see Section 6„ ) 
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Table 4 (Cont. ) 
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Table 4 (Cont. ) 
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Table 4 (Cont. ) 
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Table 4 (Cont. ) 


i) DIHEDRAL ANGLE I o 50 WAVELENGTH 5300 
AVERAGE GAIN FUNCTION U 0 E«7J 


0 

27o26 


5 

22 o 77 

- 

10 

16o26 


15 

12o44 


2 C 

12o39 


25 

12o76 


30 

1 1 o 1 9 

O 

35 

8o55 

0 

40 

6 o5 1 

<5 

45 

5 0 2A 


50 

4 o 07 




RoWoSo fluctuation 

0 

0«00 ft 


5 

o42 ft 


10 

d 28 ft 


15 

0 52 ft 


20 

o95 ft 


25 

o 28 ft 


30 

o82 ft 


35 

0 96 » 


40 

o 64 ft 


45 

o83 ft 


50 

o8l » 



j) DIHEDRAL ANGLE 2 o10 WAVELENGTH 5300 


AVERAGE GAIN FUNCTJoNU 0 E*7> 


0 

7o65 





5 

6o 13 



0 


10 

4 o4 1 


0 



13 

4 o40 


Q 



20 

5o87 



0 


23 

6o93 




0 

30 

6o92 




0 

35 

6o74 




0 

40 

6o94 




Q 

45 

4o89 




a 

50 

6o07 


RottoSo FLUCTUATION 

0 


0 

OoOO 

Q 




3 

o 16 

Q 




10 

o 1 4 





13 

o25 

0 




20 

q4 1 

a 




23 

o 17 

o 




30 

o37 

0 




33 

d?1 

0 




40 

U23 

*> 




49 

U32 

0 




50 

lo05 

0 





19 



9. RANGE CORRECTION AND PULSE SPREADING 

In order to utilize laser tracking data, the measured ranges must be extrapolated 
to the center of the satellite. The range correction that needs to be applied to the data 
is the difference between the centroid of the actual return signal and the centroid of the 
signal that would be received from a point reflector at the center of the satellite. The 
shape and centroid of the return signal vary from pulse to pulse owing to coherent inter- 
ference, satellite orientation, and position within the far-field diffraction pattern. In 
order to estimate the variation of the range correction with satellite orientation, we 
have computed the range correction for a set of 25 satellite orientations. In these 
calculations, the position of each cube comer was weighted by the total reflected 
energy from the cube corner, which is proportional to its active reflecting area. The 
average range correction for the set is 0.0733 m, with an rms deviation of 0.0007 m. 
The range correction derived in this manner is somewhat smaller than the best estimate 
obtained by weighting the position of each cube comer by the intensity of its diffraction 
pattern. 

In Table 5, the range correction, in meters, is listed in the same format and 
under the same conditions as the gain function of Table 4. For both wavelengths, the 
range correction is about 0. 075 m for the average dihedral -angle offset of 1. 5 arcsec. 

The reflected pulse from the array will be longer, on the average, than the trans- 
mitted pulse because the cube comers are not all at the same distance from the 
observer. The effect is negligible except for very short pulses, such as those from 
mode-locked lasers. One measure of the effect is the distance of the half-power point 
on the leading edge from the pulse centroid. Table 6 lists the increase of this quantity 
for the incoherent case for three different incident-pulse lengths. The values shown 
are the one-way range error that would result in a half-maximum detection system. 
Tests made with large numbers of coherent returns do not indicate that the incoherent- 
pulse spreading at the half -power point is the mean of the spreading for the coherent 
pulses. The use of Table 6 to correct range data for pulse spreading would, therefore, 
be questionable. 
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Table 5. Range correction versus velocity aberration. 
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Table 5 (Cont. ) 
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Table 5 (Cont. ) 
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Table 5 (Cont. ) 
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Table 5 (Cont 0 ) 
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Table 6. Pulse spreading. 


Pulse width 

Pulse spreading 

(nsec) 

(m) 

20.0 

0.0000 

5.0 

0.0002 

0.2 

0.0040 
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10. INFLUENCE OF OPTICAL COHERENCE 


Because of coherent interference between the reflections, the strength, centroid, 
and shape of the laser echo from the Starlette satellite vary from pulse to pulse. To 
estimate the fluctuation in measured range due to interference, we have computed sets 
of coherent pulses for three different transmitted pulse lengths (see Table 7) . The 
pulse length is given in the first column (half-power, full width), and the rms variation 
of the measured range for two types of weighting appears in columns 2 and 3. In the 
second column, all pulses are weighted equally; this procedure gives larger variations 
owing to the fact that weaker pulses tend to have greater variations in centroid. 

In the third column, each pulse is weighted by the ratio of the strength of the coherent 
pulse to the average (or incoherent) signal strength. The numbers should be considered 
only approximate since each rms variation was derived from statistical analysis of a 
set of 100 coherent returns. In general, the variations decrease as the pulse length 
decreases. 


Table 7. 

Coherent range variations. 

Pulse length 
(nsec) 

Equal 

weighting 

(m) 

Weighted by 
signal strength 
<m) 

20.0 

0.0127 

0.0072 

5.0 

0.0140 

0.0044 

0.2 

0.0068 

0.0052 


When plotted versus time, the reflected pulses from the Starlette array will have 
nearly the same shape as the transmitted pulse except for the 0. 2-nsee transmitted 
pulse length. As shown in Table 6, the reflected pulse will be significantly longer in 
this case as a result of the difference in range to individual cube comers. In addition, 
the shape of the pulse may occasionally show two peaks as a function of time, owing 
to coherent interference. 
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11. ACCURACY OF RANGE CORRECTION 


As can be seen from the reflectivity histogram of Figure 3, the spread in range of 
the reflectors contributing to the return signal is only a few centimeters. In fact, 
over 80% of the reflected energy comes from an area only 2 cm wide, and 97% is 
received from a region only 3 cm wide (see Table 3). The range uncertainty due to 
satellite geometry cannot be greater than about 1 or 1.5 cm. 

In Table 8, the range correction is given as a function of dihedral -angle offset of 
the cube corners. The numbers were obtained by averaging the values from Table 5 
over the 30- to 50-prad range of velocity aberration. 


Table 8. Range correction versus dihedral-angle offset. 


Dihedral-angle 

offset 

(arcsec) 

Wavelength 6943 

Wavelength 5300 

0 

0.0689 

0.0688 

0.5 

0.0710 

0.0705 

1.0 

0.0736 

0. 0732 

1.5 

0.0748 

0.0748 

2. 1 

0.0752 

0. 0750 


The range correction decreases as the dihedral angle decreases. None of the 
values is significantly larger than the 0. 075 m we took as the best estimate of the 
range correction. A dihedral -angle offset of 0 arcsec shows the largest difference 
(about 6 mm) from the 0. 075-m value. Since the range corrections for the larger 
dihedral-angle offsets differ by only a few millimeters, and because it is extremely 
unlikely that all the offsets are very small, the uncertainty in the computed range 
correction due to differences in beam divergence among the cube comers is probably 
less than 0. 5 cm. 
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